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Abstract 

Background: Higher seed yield is one of the objectives of jatropha breeding. However, genetic analysis of the yield 
traits has not been done in jatropha. Quantitative trait loci (QTL) mapping was conducted to identify genetic factors 
controlling growth and seed yield in jatropha, a promising biofuel crop. 

Results: A linkage map was constructed consisting of 105 SSR (simple sequence repeat) markers converged into 
1 1 linkage groups. With this map, we identified a total of 28 QTLs for 1 1 growth and seed traits using a population 
of 296 backcrossing jatropha trees. Two QTLs qTSW-5 and qTSW-7 controlling seed yield were mapped on LGs 5 
and 7 respectively, where two QTL clusters controlling yield related traits were detected harboring five and four 
QTLs respectively. These two QTL clusters were critical with pleiotropic roles in regulating plant growth and seed 
yield. Positive additive effects of the two QTLs indicated higher values for the traits conferred by the alleles from 
J. curcas, while negative additive effects of the five QTLs on LG6, controlling plant height, branch number (in the 
4th and 10th months post seed germination), female flower number and fruit number respectively, indicated 
higher values conferred by the alleles from J. integerrima. Therefore favored alleles from both the parents could be 
expected to be integrated into elite jatropha plant by further backcrossing and marker assisted selection. Efficient 
ways to improve the seed yield by applying the two QTL clusters are discussed. 

Conclusion: This study is the first report on genetic analysis of growth and seed traits with molecular markers in 
jatropha. An approach for jatropha improvement is discussed using pleiotropic QTLs, which will be likely to lead to 
initiation of molecular breeding in jatropha by integrating more markers in the QTL regions. 



Background 

Jatropha seed contains about 18-30% oil, which can be 
pressed to produce high-quality biodiesel fuel It is 
believed that jatropha can be used to alleviate the energy 
crisis, and is becoming one of the worlds key crops for 
biodiesel production [1,2]. However, jatropha must not 
directly compete with food crops due to possible food 
crisis and limited farmlands in the world, so it must be 
resistant to a high degree of aridity. Jatropha can there- 
fore be used to rehabilitate wastelands and improve the 
environment. It can also enhance the quality of rural life 



* Correspondence: genhua@tll.org.sg; chunming@tll.org.sg 
+ Equal contributors 

Molecular Population Genetics Group, Temasek Life Sciences Laboratory, 1 
Research Link, National University of Singapore, 117604 Singapore, Singapore 
Full list of author information is available at the end of the article 

(3 BioMed Central 



by providing new economic resources for marginal farm- 
lands [3]. 

Since jatropha has been domesticated, there is an im- 
mediate need to breed for superior genotypes. The 
objectives of breeding should aim at higher seed yield 
and oil content, earlier maturity, reduced plant height, 
resistance to pests and diseases, drought resistance/ 
tolerance, higher ratio of female to male flowers and 
improved fuel properties [4]. Traditional methods of 
genetic improvement of quantitative traits have relied 
mainly on phenotype and pedigree information [5], which 
are easily influenced by environmental factors. Jatropha 
seed yield is a complex trait, with difficulties in reliable 
yield prediction. In order to tackle this knowledge gap, it 
is necessary to systematically study the annual seed yield 
in operational plantation conditions along with relevant 
factors [6]. 
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As in other crops, almost all of the economically im- 
portant traits in jatropha, such as seed yield, biotic or 
abiotic stress resistance, are quantitative and determined 
by multiple genes with minor effects which are described 
as quantitative trait loci (QTL). Taking rice for example, 
a number of genes (or QTLs) for yield traits, including 
tillering [7], number of grains per panicle [8] and grain 
weight [9], have been isolated through map based clon- 
ing. The genes regulating yield traits and developmental 
processes can function at various stages, in different 
pathways and through diverse mechanisms in rice [10]. 
A genomic region on chromosome 4 was detected with 
multiple effects on increased flag leaf width and length, 
and panicle number and length in rice [11]. Clustered 
QTLs were recently reported for source leaf size and 
yield traits in rice [12], and improving rice yield and 
quality by QTL pyramiding was carried out [13]. Marker 
assisted breeding has been applied by taking advantages 
of the useful information on these genes or QTLs affect- 
ing agronomic traits of importance. 

In contrast, jatropha, which was still considered wild 
in 2010 [6], had not yet undergone a careful breeding 
program with systematic selection and improvement of 
suitable germplasm. Recently, we have established a first 
generation genetic linkage map using 506 microsatellite 
and SNP (Single Nucleotide Polymorphism) markers 
covering 11 linkage groups [14], and conducted a whole 
genome scan for QTL and eQTL affecting seed oil traits 
[15]. In addition, we have isolated and identified miR- 
NAs and targets in jatropha [16]. However, the genetic 
bases of jatropha growth and seed yield have not been 
studied in jatropha breeding. 

In this paper, we describe the genetic bases of seed 
yield through QTL mapping, which is one of the most 



important agronomic traits, together with plant height, 
stem diameter, branch number, female flower number 
and fruit number. Moreover we have analyzed the pleio- 
tropic effects and interaction of the QTLs, and provided 
an approach for possible modulation of the QTLs to im- 
prove growth and seed character in jatropha. 

Results 

Trait analysis 

Growth and seed traits were measured in a QTL map- 
ping population, and the frequency distributions of all 
traits in the progeny showed a continuous distribution. 
The distribution of phenotypic values showed bi- 
directional transgressive segregation (Table 1), revealing 
complex genetic bases of these traits. While seed yield in 
/. curcas was higher than that in /. integerrima, branch 
number in /. integerrima is significantly higher than that 
in /. curcas. The data implied that /. integerrima germ- 
plasm could be applied for hybrid breeding to improve 
agronomic traits, such as branch number in the fourth 
and tenth months, and female flower number. 

Correlation analysis among these traits was performed 
(Table 2), and total seed weight showed a significant cor- 
relation with total branch number, female flower number 
and fruit number, with coefficients 0.364, 0.294 and 
0.308, respectively. Therefore, these agronomic traits 
were suggested to be key factors for seed yields. 

QTL mapping 

The linkage map consisting of 105 DNA markers and 
covering 643.8 cM of the genome, converged into 11 
LGs (linkage groups) corresponding to 11 chromosome 
pairs in jatropha. The average distance between markers 



Table 1 Descriptive statistics on phenotype data of QTL mapping population and parents (J. curcas PZMD16, 
J. integerrima S001 and F1 CI7041) 



N 


Trait 


Acronym 


Mean 


SD 


Min 


Max 


PZMD16 


CI7041 


Growth traits 


1 


Height in the 4th month 


H4M 


70.5 


24.6 


15.0 


134.0 


62.16 


70 


2 


Height in the 10th month 


H10M 


122.9 


37.2 


33.0 


272.0 


166 


180 


3 


Diameter in the 4th month 


D4M 


1.9 


0.5 


0.6 


3.2 


1.45 


1.2 


4 


Diameter in the 10th month 


D10M 


4.3 


1.1 


1.1 


7.0 


6.2 


5 


5 


Branch number in the 4th month 


BN4M 


1.8 


2.2 


0.0 


14.0 


10 


12 


6 


Branch number in the 10th month 


BN10M 


6.4 


4.0 


1.0 


19.0 


15 


17 


7 


Total branch number 


TBN 


4.4 


2.1 


1.0 


14.0 


3.5 


3.33 


8 


New branch number per branch 


BNPB 


2.9 


1.1 


1.0 


8.0 


1.75 


1.67 


Flower, fruit and seed yield 


9 


Female flower number 


FFN 


4.5 


2.8 


0.0 


15.0 


4.75 


9.25 


10 


Fruit number 


FRUITNO 


9.1 


10.5 


0.4 


62.0 


35 


No fruit 
(Hybrid F1) 


11 


Total seed weight in 2010 


TSW 


34.1 


60.1 


0.5 


541.4 


360.0 


90.0 
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Table 2 Correlation coefficients and significance of correlations among growth and yield traits in a QTL mapping 
population 



H4m HlOm 


D4m 




DIOm 


BN4m 


BNIOm 


TBN 




BNPB 


FFN 


FruitNo 


H10m 0.371*** 


D4m 0.652*** 0.374*** 


D10m 0.171*** 0.675*** 


0.494*** 




















BN4m 0.170*** 0.064 


0.229*** 




0.177*** 
















BN10m 0.199*** 0.189*** 


0.175*** 




0.401*** 


0.415*** 














TBN -0.155* 0.041 


-0.092 




0.154* 


0.067 


0.006 












BNPB -0.027 -0.015 


-0.005 




0.094 


0.141 


0.072 


0.450* 










FFN -0.448*** 0.138* 


-0.207*** 




0.258*** 


-0.042 


-0.061 


0.265* 




0.097 






FruitNo 0.026 0.220*** 


0.077 




0.256*** 


0.198** 


0.327*** 


0.097 




0.038 


0.215** 




TSW -0.137* 0.093 


0.004 




0.262*** 


0.156* 


0.164** 


0.364*** 


0.048 


0.294*** 


0.308*** 


P values are as follows: * P < 0.1 0, ** P < 0.05, *** P < 0.01 . 




















Table 3 QTLs for growth traits, seed characters 
















Trait 


QTL a 


Linkage Group 


IVId 1 Ktri 


Prtcitirtn ^ 
rUbl UUI 1 


cM LOD Peak R 2 c (%) 


Additive Effects d 


Growth H4m 


qH4m-3 


3 




In iint1 7Q 
JLUII IL I / y 


JJ.O 




3.50 


5 




-17.51 




qH4m-7 


7 




Jd Ll O 1 U 


71 Q 
/ I .y 




2.97 


4.4 




-10.31 


H10m 


qH10m-3 


3 




Jd Ll I UJ^- 


1 1 Q 
I I .y 




3.55 


6.3 




21.53 




qH10m-5 


5 




Jd Ll y^J 


42 1 




5.81 


8.5 




21.99 




qH10m-6 


6 




In iint3 1 1 
JLUII ILJ 1 Z 


Z J.J 




2.56 


4 




-15.80 




qH10m-9 


9 




Jd Ll OJy 


J.U 




3.09 


4.5 




16.30 




qH10m-10 


10 




In iintflR1 
JLUII ILUo 1 


ZU.D 




4.63 


6.7 




-19.25 


D4m 


qD4m-5 


5 




wlcl 


J I .J 




3.88 


5.6 




0.26 




qD4m-9 


9 




Jd Ll KJyO 


24 6 




2.70 


5.2 




0.25 




qD4m-11 


11 




Jd Ll DO^L 


1 ZL 3 




2.98 


4.6 




0.25 


D10m 


qDWm-30 


3 




Jd Ll I UJ^L 


1 1 Q 
I I .y 




8.11 


12.7 




0.96 




qDWm-3b 


3 




ln lintflAR 
JLUII ILU^-O 


c;q Q 




3.85 


4.9 




-0.68 




qD10m-5 


5 




Jd Ll /HO 


3R 3 
Jo. J 




15.03 


21.1 




1.11 


BN4m 


qBN4m-6 


6 




In lintO^ft 
JLUII ILUJO 


OHO 




4.04 


6.9 




-1.26 


BN10m 


qBN10m-l 


1 




JdLI / ZZ 


^ZL 1 

J-H-. I 




3.44 


5.6 




2.04 




qBN10m-4 


4 




JdLIOJ^f 


zli n 




3.51 


5.5 




-2.19 




qBN10m-6 


6 




Jcuintl 1 1 


93.4 




3.58 


5.9 




-2.28 


TBN 


qTBN-7 


7 




Jatr610 


68.9 




3.40 


8.4 




1.30 


BNPB 


qBNPB-5 


5 




Jatr739 


46.2 




2.98 


7.9 




0.62 


Flower, Fruit and seed yield FFN 


qFFN-2 


2 




Jatr691 


20.7 




2.00 


4.5 




1.20 




qFFN-5 


5 




Olel 


30.3 




3.38 


7 




1.56 




qFFN-6 


6 




Jatr301 


15.8 




6.41 


13.6 




-2.21 




qFFN-7 


7 




Jcuintl 51 


2.0 




4.16 


9.6 




1.97 


FruitNo 


qFruitNo-1 


1 




Jatr749 


42.4 




2.47 


5.5 




5.54 




qFruitNo-6 


6 




Jatr839 


28.3 




4.97 


12.4 




-7.61 




qFruitNo-7 


7 




Jatr866 


64.8 




3.15 


7.3 




5.65 


TSW 


qTSW-5 


5 




Jcuint002 


33.7 




2.24 


5.2 




29.33 




qTSW-7 


7 




Jatr866 


66.8 




2.70 


4.9 




26.31 



a QTL: starting with "q," followed by an abbreviation of the trait name, the name of the linkage group, and the number of QTL affecting the trait on the linkage 
group. b Position from the first marker on each linkage group. c Coefficient of determination or the percentage of variance explained (PVE) by the detected QTL. d 
Estimated phenotypic effect of substituting J. integerrima alleles with J. curcas alleles at QTL. 
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was 6.6 cM. Most of the LGs were consistent with those 
described previously [14]. 

QTL analyses were performed on the means of growth 
traits, branch number, female flower and fruit number, and 
seed yield (Table 3; Figure 1). We have detected 28 QTLs 
for all traits examined with LOD threshold 2.0 to 2.5 deter- 
mined by permutations. Individual QTLs were detected 
with percentage of variation explained (PVE or R 2 ) 3 to 
21.16%, and four of them had PVE exceeding 10%. 

QTLs with positive and negative additive effects were 
identified, with a positive effect implying a higher value 



for the trait conferred by the allele from /. curcas, and 
negative from /. integerrima (Table 3). 

QTLs for growth traits 

Sixteen QTLs were identified and dispersed among all 
the linkage groups except LGs 2 and 8. Four QTLs over- 
lapping on the lower part of LG5, namely qHlOm-5, 
qD4m-S, qDlOm-5 and qTBN-5, were detected under- 
lying plant height in the 10th month, stem diameter in 
the 4th and 10th months, and total branch number, re- 
spectively (Figure 1). Additive effects of these QTLs 



LG1 



LG6 



16.9 ^sj/^ 



19.2 — 
20.6- 
21.3 V 



/~ Jcuint098 
r Jatr301 

- Jatr773 
^Jatr1013 

- Jcuint312 
~- Jcuint311 




LG2 




LG3 



0.0 — P| — Jatr856 

I 



8.9 — H — Jatr1054 
14.9 M KASI 



- Jcuint105 

- Jcuint143 
HCurcinl Curcin2 




LG7 



0.0 — f— ) — Jcuint151 



LG8 



T 




43.3- 
50.2- 



LG4 



LG5 




Jcuint094 
Jatr735 

Jcuint354 Jatr945 
Jatr739 



i! 



II 



It 



LG9 



LG10 



LG11 



-Jcuint012 0.0- 




Jcuint131 Jcuint081 
Jcuint067 
Jcuint066 



- Jcuint044 

- Jcuint020 
-KASII 



- Jatr684 

- Jatr306 



Figure 1 Summary of QTL locations detected. QTL represented by bars are shown on the left of the linkage groups, close to their 
corresponding markers. The lengths of the bars are proportional to the confidence intervals of the corresponding QTL in which the inner line 
indicates position of maximum LOD score. The confidence intervals of QTL are showed in green color, and two QTL clusters are highlighted in 
red. 
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were positive, indicating that the alleles from /. curcas 
increased these trait values. 

Conversely, two QTLs, namely qBN4m-6 and 
qBN10rn-6, were detected on the lower part of LG6 con- 
trolling branch number with negative additive values, in- 
dicating Jatropha integerrima allele increased branch 
number. 

QTLs for female flower and fruit number 

Six QTLs were identified and dispersed on LGs 1, 5, 6 
and 7, with two QTLs, namely qFFN-6 and FruitNo-6, 
being located on the same region of LG6, controlling fe- 
male flower number and fruit number respectively. The 
PVE of these two QTLs were higher than 10%, indicating 
their significant effects on the two important yield trait 
components. 

QTLs for seed traits 

On LGs 5 and 7, two QTLs of qTSW-S and qTSW-7 
were detected controlling total seed weight, which is one 
of the most economically important traits. Interestingly, 
QTLs underlying yield related traits were clustered at 
these two QTLs. At qWT-5, four QTLs underlying plant 
height, stem diameter, branch number and female flower 
number were detected. Near qTSW-7, three QTLs of 
qH4m-7, qTBN-7 and qFruitNo-7 were detected, con- 
trolling plant height, total branch number and fruit 
number respectively. 

It was noteworthy that two QTL clusters were detected 
on LGs 5 and 7, respectively. Five QTLs were detected 
on the lower part of LG5 (Figure 2A), and four QTL clus- 
ters were detected on lower part of LG7 (Figure 2B). 

Favored alleles originated from two parents 

Two QTL clusters were detected consisting of five and 
four QTLs, controlling total seed weight, plant height, 
stem diameter, female flower number and fruit number. 
The positive additive effects indicated higher values for 
the traits conferred by the allele from /. curcas. Mean- 
while five QTLs on LG6, namely qH4m-6, qBN4m-6, 
qBN10m-6 qFFN-6 and qFruitNo-6, controlling plant 
height, branch number (in 4th and 10th months post 
seed germination), female flower number and fruit num- 
ber respectively, were detected with negative additive 
effects indicating higher values conferred by /. integer- 
rima (Table 3). 

Major effects of qTSW-5 and qTSW-7 

A two-way analysis of variance (ANOVA) was carried 
out to assess genetic effects and interactions of the two 
QTLs of qTSW-S and qTSW-7 controlling total seed 
weight. The values of different genotypes are shown in 
Figure 3. Total seed weight was significantly increased in 
the presence of these two QTLs. When qTSW-S 



presented, total seed weight was improved from 
16.66 ±7.26 to 42.00 ±5.06 g, and qTSW-7, from 
15.97 ± 6.36 to 42.69 ± 6.16 g (Figure 3A). 

Interestingly, we found that the two QTLs for seed 
yield overlapped with other QTLs for other agronomic 
traits than seed yield itself. ANOVA showed that the 
QTL qTSW-S for seed yield affected significantly plant 
height, stem diameter, new branch number per branch 
and female flower number, while qTSW-7 affected plant 
height, total branch number and fruit number 
(Figure 3B). 

Effect of pyramiding qTSW-5 and qTSW-7 

The interaction between marker effects for qTSW-S and 
qTSW-7 was non- significant with a relatively low P value 
(0.14) (Table 4), while the marker effects for qTSW-S 
and qTSW-7 were non-additive (Figure 4). This could be 
caused by the lack of power in the ANOVA due to an 
unequal distribution of genotypic classes (Figure 3). 

Despite the non-significance of the interaction of the 
two QTLs, total seed weight was significantly increased 
in the presence of the two QTLs. Lines carrying both 
QTLs produced an average 61.93 ±7.31 g of seeds, 
nearly three times as much as any other genotype com- 
binations (Figure 4). Therefore, although total seed 
weight could be improved by introducing the two QTLs, 
there would be advantages to be gained by pyramiding 
the two QTLs. 

Discussion 

Pleiotropic QTLs for growth and seed yield 

Improvement of yield potential is the most important 
goal of most breeding programs worldwide. However, 
yield is a complex trait controlled by many genes with 
major or minor effect [17]. Till date, selective breeding 
for high yield remains the most challenging task. Our 
results show that QTL clusters may have multiple effects 
on yield and yield-related traits, and we have detected 
two QTL clusters controlling multiple traits on the same 
regions of LGs 5 and 7. The two QTLs for seed yield 
clustered with those for plant height, stem diameter, 
branch number, female flower number or fruit number, 
revealing that these are two critical regions for jatropha 
growth and seed yield. Similarly, a major QTL, Ghd8, 
plays pleiotropic roles in regulating grain productivity, 
plant height, and heading date in rice [18]. These yield- 
related traits, such as height, diameter, branch number, 
female flower number, fruit number and so on, are less 
complex than total seed yield per se, yet highly corre- 
lated with total seed yield, hence it could be very useful 
for selection at different breeding stages for evaluating 
respective agronomic traits [19]. 

Some of the pleiotropy could be a consequence of cor- 
relations among traits such as female flower number, 
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A 

16C* LC0 



12 8- - 



9.6 - • 



6 .4 - - 



3 2 - - 



0 0 



H4m 

H1 Om 

D4m 
D1 Om 
BN4m 
BN1 Om 
TBN 
BNPB 
FFN 
FruitNo 
TSW 



/ 

\ 

: 

u 

if 



/ 



/ 



/ 



/ 



A 

/ X 




10 



is 



20 25 
Chromdsom©-5 (LG5) 



30 



36 



40 



Jcuint002 



46 



51 



B 

5 0 



4 0 - • 



3.0 - • 



2 0 - * 



10- 



00- 



LODO 




\ 




15 



22 29 36 

Chromosome-7 (LG7) 



58 J 66 

Jatr866 



73 



Figure 2 QTL clusters on LGs 5 and 7. QTL scans of growth on linkage maps. Horizontal line indicates 5% LOD significance thresholds (2.0) 
based on permutation. A: LG5; B: LG7. 
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A 



60.00 



50.00 
% 40.00 

30.00 
J 20.00 
£ 10.00 
0.00 



42.00 



42.69 



16.66 



15.97 



Aa AA 
qTSW5 



Bb BB 
qTSW7 



B 



132.00 
130.00 
128.00 
126.00 
124.00 
122.00 
120.00 
118.00 
116.00 
114.00 
112.00 
110.00 
8.00 
6.00 
4 00 
2.00 
0 00 



114.57 



128.79 



4 69 4 .96 

3.73 _ . ,nr 3.91 



» <■* n n ~ n 1*1 

n □ n II n n 



80.00 73.07 
70.00 
60.00 
50.00 
40.00 
30.00 
20.00 
10.00 
0.00 



66.97 



3.74 4.92 6 - 27 



11.61 



Aa | AA 


Aa AA 


Aa | AA 


Aa | AA 


Aa | AA 


H10m (cm) 


D4m (cm) 


D10m (cm) 


BNPB 


FFN 








(number) 


(number) 



qTSW-5 



Bb 



BB 



Bb 



BB 



Bb 



n 

BB 



H4m (cm) TBN (number) 



qTSW-7 



FruitNo 
(number) 



Figure 3 Total seed weight (g) and related traits of plants with different genotypes of qTSW-5 (AA, Aa) and qTSW-7 (BB, Bb); N denotes 
sample number of each genotypic classes; Error bars denote SEs (Standard Errors). A: Significant major effects of the two QTLs on seed 
yield; B: The two QTLs with pleiotropic roles in regulating plant growth and seed yield. Significant at P < 0.01 of Bonferroni test. 



fruit number and seed yield, or between diameter and 
height, which were highly correlated at the phenotypic 
level Some of the co-occurrence of QTLs could be a 
natural consequence of allometry, as has been suggested 
in the study on poplar [20]. In spite of that, the allomet- 
ric relationship may not exist in alternative genetic back- 
grounds that have been exposed to different selection 
pressures [20]. 



The pleiotropic QTLs could be explained in different 
ways. Chromosomal regions were associated with more 
than two traits indicating either linkage or pleiotropic 
effects. There could be certain genes coexisting in these 
QTLs or a certain gene with pleiotropic effects on j atro- 
phia growth and seed development. Here, the QTLs we 
reported were still distant to the flanking markers with 
linkage distance, therefore, it will be meaningful to 
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Table 4 ANOVA of seed yield in the QTL mapping 
population based on genotypes of the marker loci that 
are most closely linked to the QTLs 



Effect 


d.f. 


MS 


F 


P 


1 (Jcuint002, qTSW-5) 


1 


24943.135 


8.190 


0.005** 


2 (Jatr866, qTSW-7) 


1 


27739.644 


9.110 


0.003** 


1 X 2 


1 


6703.480 


2.200 


0.140 



** Denotes significant at P = 0.01. 

conduct fine mapping of these QTLs, isolate the target 
genes, and understand whether linkage or pleiotropic 
effects are responsible. Fine mapped QTL will speed up 
genetic improvement through marker assisted selection 
(MAS) by applying the closely linked markers [21]. As 
molecular markers are still limited in jatropha, we are 
constructing a second generation linkage map of jatro- 
pha with a high-resolution of SSR or SNP markers, 
which will lay a solid foundation for a variety of future 
genetic and genomic studies, including QTL fine map- 
ping and marker assisted breeding. 

Towards molecular breeding by transferring favored 
alleles from the two parents 

Plant growth and seed traits in jatropha are controlled 
by multiple gene complexes. Genetic markers have made 
it possible to detect QTLs that are significantly asso- 
ciated with traits, making selection more effective [21]. 
Genetic response can be further improved by inclusion 
of the QTLs in marker assisted breeding, which makes 
use of phenotypic, genotypic and pedigree data [22]. 
DNA markers have enormous potential to improve via 
MAS the efficiency and precision of conventional plant 
breeding [23], including jatropha germplasm enhance- 
ment and genetic improvement. The exploitation of the 
advantages of MAS relative to conventional breeding 
could have a great impact on crop improvement. We 
have identified markers linked to some major QTLs and 
genes by constructing a backcross population between 
/. curcas and /. integer imma. Further interspecific 
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Figure 4 Effects of pyramiding the two QTLs of qTSW-5 
(genotypes of AA, Aa) and qTSW-7 (BB, Bb) on seed yield. Error 
bars denote SEs. 
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AA 


AA 
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N(93) 
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introgressions could be expected to apply the allelic 
sources for trait improvement. 

The QTLs of qTSW-5 and qTSW-7 controlling seed 
yield were detected in two QTL clusters on LGs 5 and 7 
respectively. The positive additive effects indicated 
higher values for the traits conferred by the alleles from 
/. curcas. Meanwhile the five QTLs on LG6, controlling 
plant height, branch number (in the 4th and 10th 
months post seed germination), female flower number 
and fruit number, were detected with negative additive 
effects, indicating higher values conferred by /. integer- 
rima. Therefore, it will be feasible to transfer favored 
alleles from both the parents to elite jatropha varieties as 
recurrent lines. Hence, the QTL mapping population 
will be very useful in transferring favored alleles from 
both the parents by further backcrossing and marker 
assisted breeding. 

QTL pyramiding 

Accumulating major genes for seed yield in an elite 
genotype by conventional breeding is laborious and 
time-consuming [24]. Gene pyramiding is difficult using 
conventional phenotyping methods, due to the epistatic 
effects of genes [25]. However, an advantage of pyramid- 
ing the two QTLs governing seed yield was observed in 
our study, and the identification of markers linked with 
each QTL allows for the identification of plants carrying 
one or both QTLs. In rice, rice yield and quality were 
improved by QTL pyramiding. The pyramid line 
(qHD8 + GS3) had higher yield potential, longer grains, 
and a more suitable heading date [13]. Our results reveal 
that detailed information regarding the pyramiding effect 
is very important for efficiency of marker assisted pyra- 
miding of different alleles at target QTL in jatropha 
breeding. The markers described here may serve as use- 
ful tools for gene pyramiding with the two QTLs. 
Advanced lines with a good genetic background and 
high seed yield genes combinations will be expected to 
have great practical breeding value. Nevertheless, two 
issues need to be further addressed, one is whether 
effects of the QTLs are family specific; another is 
whether these QTLs will be robust enough to apply to 
other environments. Therefore the markers need to be 
further tested in different families and environments. 

Conclusions 

We identified a total of 28 QTLs underlying the growth 
and seed yield traits in jatropha. This study represented 
the first investigation on plant growth and seed yield 
through QTL mapping in jatropha. An approach was 
discussed for jatropha improvement using pleiotropic 
QTLs, which could be likely to lead to initiation of MAS 
by integrating more markers in the critical regions of the 
two QTL clusters. 
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Methods 

Plant material and growth conditions 

/. curcas PZMD16 was crossed to /. integerrima S001 
and hybrids Fl lines were generated. Then a BC1F1 
population was constructed consisting of 296 individuals 
derived from the backcross between PZMD16 used as 
recurrent parent line and an Fl line named as CI704L 
The population and parental lines were planted under 
standard growth conditions with spacing 2 m X 2 m in 
2008 in experimental field of Lim Chu Kang farm, Singa- 
pore. Due to jatrophas perennial life cycle, we generated 
backcross populations in 2008 and observed phenotypes 
from 2009. 

We used the fertilizers 15-15-15(N-15%, P205-15%, 
K2(M5%) and 13-13-21(N-13%, P205-13%, K20-2196) 
to promote growth of flowers and fruits. The fertilizer 
application was three times per year. In 2008, we applied 
100 g of each fertilizer per plant each time. And we ap- 
plied 200 g of each fertilizer per plant each time from 
2009. The insecticides Rogor L-40 (Dimethoate 38% w/ 
w) or Alcalineum (Mineral 80% w/w) were used in the 
farm every two months with the concentration of 3-5 ml 
/10 L. The fungicides were Mancozide WP (Mancozeb 
80% w/w) or SAPROL(Triforine 17.8% w/w), applied 
once every two months with the concentration of 22 g/ 
10 L water or 12.5 ml/10 L water, respectively. 

DNA markers and genotyping 

We extracted total DNA from leaves using the DNeasy 
plant mini kit (QIAGEN, Germany). We selected one 
hundred and five markers almost evenly covering the 11 
LGs from a first-generation linkage map of jatropha 
[14]. One primer of the selected markers was labeled 
with FAM or HEX fluorescent dyes at the 5' end. The 
PCR for microsatellite amplifications on PTC- 100 PCR 
machines (MJ Research, CA, USA) was done using the 
program as follows: 94 °C for 2 min followed by 37 
cycles of 94 °C for 30 s, 55 °C for 30 s and 72 °C for 
45 s, then a final step of 72 °C for 5 min. Each PCR reac- 
tion consisted of lx PCR buffer (Finnzymes, Espoo, Fin- 
land) with 1.5 mM MgCl 2 , 200 nM of each PCR primer, 
50 uM of each dNTP, 10 ng genomic DNA and one unit 
of DNA-polymerase (Finnzymes, Espoo, Finland). Pro- 
ducts were analyzed using a DNA sequencer ABI3730xl 
(Applied Biosystems, CA, USA), and fragment sizes were 
determined against the size standard ROX-500 (Applied 
Biosystems, CA, USA) with software GeneMapper V4.1 
(Applied Biosystems, CA, USA) as described previously 
[26]. 

Agronomic trait measurement and data collection 

The growth traits were observed in the 4th and 10th 
months post seed germination because most of the 
plants started flowering in the 4th month and completed 



growth in the 10th month. We decided to observe the 
growth traits in the two critical times. Phenotypic data 
were collected from the QTL mapping population as 
follows: 

Growth traits included height, diameter, branch num- 
ber in the 4th and 10th months, total branch number 
and new branch number per pruned branch. Pruning is 
an agronomic treatment to produce more branches. 

Flower and fruit traits included female flower number 
and fruit number, which were observed in three flower 
and fruit clusters respectively. The average numbers of 
the three replications were used for further QTL 
analysis. 

In 2010, we harvested the fruits of each plant with sep- 
arate nylon mesh bags twice a month and dried them in 
the sun. Then we hulled them to get the seeds. After 
recording the seeds' number and dry weight, we stored 
them inside a 4 °C cooling room. All the seeds harvested 
in 2010 were collected for evaluating total seed weight 
of one year. Seeds were harvested at maturity after the 
color of the fruits had changed from green to yellow- 
brown. 

Statistical analysis 

Plant height, diameter and branch number traits were 
collected at different stages in the backcross population 
consisting of 296 individuals, and traits of flower and 
fruit were collected with three replications. Pearson 
phenotypic correlations were calculated among all the 
traits using SAS PROC CORR [27]. 

Linkage map was constructed using the software 
CRIMAP 3.0 [28] with the genotyping data of 105 mar- 
kers in the QTL mapping population. Kosambi function 
was used to calculate all multipoint distances. Graphical 
visualization of the linkage groups was completed with 
MapChart 2.2 software [29]. QTL analysis was carried 
out using QTL Cartographer version 2.5 [30]. Composite 
interval mapping (CIM) was used for mapping QTLs 
and estimating their effects. The forward regression 
method was used to scan the genome at 2-cM intervals. 
The log of the odds (LOD) score was determined for de- 
claring a significant QTL by permutation test analyses 
(1,000 permutations, 5% overall error level). 

The position and confidence interval of QTL were 
determined as described previously [15]. Briefly, the 
maximum LOD score was taken as the position of the 
QTL, and the region in the LOD score within 1 LOD 
unit of maximum was taken as the confidence interval. 
Additive effects of the detected QTL were estimated as 
the mean effects of replacing hybrid (CI7041) s alleles at 
the locus of interest by /. curcas (PZMD16) alleles. Thus, 
for a QTL to have a positive effect, the /. curcas alleles 
must increase the trait value. The contribution of each 
identified QTL to total phenotypic variance (R 2 ) or 
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percentage of variance explained by the QTL (PVE) was 
estimated by variance component analysis. QTL nomen- 
clature was adapted as follows: starting with "q," fol- 
lowed by an abbreviation of the trait name, the name of 
the linkage group and the number of QTL affecting the 
trait on the linkage group. 

Using marker genotypes as the groups, analysis of vari- 
ance (ANOVA) was performed with the general linear 
model (GLM) procedure of SAS [27]. 
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